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INTRODUCTION 


This manual presents a general description of the STARS -2B digital V 
computer program. FORTRAN IV is used exclusively in writing the various 
subroutines. The execution of this program requires the use of eleven 
temporary storage units. 

The program was initially written and debugged on the IBM 370-165 computer 
and then converted to the UNIVAC 1108 computer, where it utilizes 
approximately 46,000 words of core. Only basic FORTRAN Library routines 
are required by the program, these being: sine, cosine, absolute value, 
square root, min 0, max 0, amax 1, anain 1, and random. 

For ease and speed in usage, the Table of Contents on the following page 
has also been laid out to present the call sequence of the program. 
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SUBROUTINE MAIN 


MAIN is the control link for the entire program. Sizing values are read 
into the program, as well as information for eigenvalue or determinant 
accuracy calculations, and the material property tables. Calls are made 
to subroutines RIEMAN and SEGMAT once for each segment in a region; then 
subroutine REGMAT is called. This procedure is executed once for every 
region in the structure. Finally a call to . subroutine STRMAT is made. 

If the clues are set for the program to calculate the prestress state, _ 

• i .V 

subroutines INITAL and LEBEGE are called. Otherwise the control link 
loops back and resets the clues for the calculation of nonlinear stiffiiess 
matrices . 

In the calculation loop for the nonlinear stiffness matrices the program 
proceeds as above to STRMAT. Next EIGVAL or DETERM are called to calculate 
eigenvalues or evaluate a determinant as per user input. If these calculations • 
satisfy convergence criteria, subroutines BCVECT, INITAL, and LEBEGE are 
called. If eonvegence criteria are not satisfied the program will loop until 
satisfaction is obtained. 

i 

In a multi-harmonic eigenvalue search the whole procedure loops on the 
harmonic number. 

There are also several counters in this control link. These are defined as 

I 

follows: . 

» 

NSC - Counts the calls to subroutines RIEMAN and SEGMAT, from 
1 to the , number of segments within a region. 

NRC - Counts the calls to subroutine REGMAT, from 1 to the 
number of regions in the structure. 

ICYC - Cycle counter (counts harmonics). 

NPASS - Pass counter (counts passes necessary to converge in a 
single harmonic). 

The block data and overlay listings are included in this section. 



READ: 
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SUBROUTINE RIEMAN 


This subroutine link assembles the data tables for use in the integration 
procedure. The program has the capability of handling 2 loading conditions 
(dead and live). Temperature load can be included in either one, but 
not both. 

The subprogram link, RIEMAN, utilizes the subroutines SETUP, ROBOT, DIF1 
or DIEF2, to integrate the differential equations of each segment 
independently, under arbitrary load conditions. The results of the 
integrations of each segment are stored in the YCORR array in RIEMAN, and 
represent the stiffness and -deflection coefficients of each segment. 
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SUBROUTINE SETUP 


SETUP is a double entry subroutine called from RIEMAN. It is a mixed precision, 
numerical integration routine, with automatic selection of a variable integration 
step size, which utilizes fifth order Runge-Kutta equations to obtain the solution 
for first order differential equations. 

SUBROUTINE MAGIC 

MAGIC is an alternate entry point to subroutine SETUP. 
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RUNGE-KUTTA GENERAL FLOW CHART 


ENTRY SETUP 


INITIALIZE MAGIC 
CALL SETUP (K = 0) 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. 1.) 

RO 

r o 

R1 • 

' r l 

R1DOT 

r l ,4- 

CS 

cos <j> 

SN ' " ' ' - 

sin <f> 

A 

a 

C 

c 

XN 

n 

F2 

f 2 

F3 

f 3 

TAN; TN ■ 

t an <j> 

SEC 

sec <j> 

Til 

T. . 
11 

TIK 

T. 

1C 

TOK 

T 

°c 

TOO 

T 

oo 

TEFREE 

T 

HI 

h. 

l 

HO 

h 

, o 

T 

t 

TI 

t . 
i 

TO 

t 

o 

SNSQ 

sin <j> 

. CSSQ 

2 

COS. <p 

CN 

cos 4) sin 4> 

X1CS 

l/cos <p 

X1SN 

l/sin 4> 

R2 

r 2 

BETA 

P 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. 1 ) 

X1ROSN 

l/r Q sin <J> 

X1ROCS 

l/r Q cos <j) 

CSX1RO 

cos 4>/r Q 

CSX1R1 

cos 4>/r^ 

CSX1R2 

cos <p/r 2 

SNX1RO 

sin <p/r Q 

SNX1R1 

sin (f>/r^ 

X1R1 

l/r ! . 

X1R2 

l/r 2 

X1R1SQ 

i/ri 2 

X1ROSQ 

OJ 

o 
1 — 1 




29 





ROBOT 


ENTRY 
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SUBROUTINES DIFl AND DIFF2 


These subroutines are called in RIEMAN as necessary. DIFl contains 
the differential equations for the THIC and ST clues, while DIFF2 contains 
the differential equations for the RWA and ISG clues. Geometry clues, 
trigonometric values, and predicted values of the differential equation 
variables are passed via label common area, EQDAZN, to subroutines DIFl 
or DIFF2. The coefficients for nonlinear and load terms, XI, X2 and K, 
are identified depending upon the pass number and convergence criterion. 

The specific derivative equations and auxiliary equations are contained in 
these subroutines. The values of each derivative equation, YDOT, and each 
auxiliary equation, YA — are returned to RIEMAN via label common EQUAZN. 

i 

A special equation counter, I, is used in these subroutines, which counts 
in increments of eight. The first eight values of I, 1 through 57 ( in 
increments of eight), correspond to the eight sets of initial conditions 
required to compute the segment stiffness matrices in subroutine SEGMAT. 

The subsequent values of I, 65 and 73 maximum (again in increments of 
eight) correspond to the computation of each set of eight equations for 
each loading condition (two conditions corresponding to dead and live loads 
respectively). 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. 

i ! 

XN 

n 

dT 

YDOT (I ) ' 

T „ 

<t> 

(J)0 ,<j> 

ds 



dN , 

YDOT (I + 1) 

N 

i 


ds 



dJ, 

YDOT (1+2) 

J, 

<P 

<P ,<t> 

ds 



dM 

YDOT (1+3) 

M , 

4> 

•• 

<P »<P 

ds 

YDOT (1+4) 

U, , 

dU 


ds 

YDOT' (1 + 5) 

V, A 

dV 

’<f> . 

ds 

YDOT '('I +’ 6) ' ' ' 

' ' w. 

dW 

*4> 

ds 

■ , 


d % 

YDOT (I + T) 


e »4> 

ds 

YPRED (I) 

T 

4>e 


YPRED (I + 1) 

N <! > 


YPRED (1+2) 

J <{> 


YPRED (1+3) 

M * 


YPRED (1+4) 

U 


YPRED (I + 5) 

V 


YPRED (I + 6) 

w 


YPRED (1+7) 

n e 


YAMPT 

CD 

s 


YANTH 

N e 


YAMTH 

! M e 
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FORTRAN CODE 

ERGIREERIRG SYMBOLS (REF. i ) 

R2SQ 

2 

r 2 

ROSQ 

2 ■ 
r o 

X1R0' 

O 

U 

i — 1 

S 

■ s ; 

XK12 

K 12 

i 

XK21 

1 — I 
CM 

XD12 

D 12 

XD21 

D 21 

XC11 

C 11 

XC22 

C 22 

XNSQ 

2 

n 


Ron-Linear Redefinitions (Ref. 1 ) 


YDOT (1+2) 

YPRED (1+2) 

YAJPH 

XRL 

XNPHI 


*J 


^ 9 *1 


»J, 


a 

N. 


d*J, 


ds 
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DIF 1 
ENTRY 


1 


IDENTIFY K, XI, 
X2 FOR EACH 
ANALYSIS PASS 


IF KGEOM = 
IGEOM = 1 

, 2, 5, 6, OR 7 

' 


IF KGEOM = 

3; IGEOM = 2 

1 

f 

IF KGEOM = 

4; IGEOM = 3 
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SUBROUTINE SEGMAT 


The results of the subroutine link, RIEMAN, are passed through the label 
common area, IZCORR, to this subroutine. SEGMAT places the elements of 
the YCORR array into several double-subscripted arrays, forms some ■ 
coordinate transformation arrays, and calls subroutine SREVN2 for double 
precision matrix inversion. 

As a result of appropriate matrix operations this subroutine produces a 
segment stiffness matrix, the XKS array, and a segment load matrix, the 
XLS array, for each segment. In passes other than the first pass, first 
cycle, the XLS array is not calculated. SEGMAT also orients each segment 
into the global coordinate system of the structure as a result of the 
matrix operations. 

Subroutine SREVN2 

SREVN2 is a subroutine called by SEGMAT to invert a real, double -precis ion, 
in-core matrix utilizing Gauss-Jordan elimination with partial pivoting. 
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FORTRAN CODE 

SNI 
SNJ 
CSI ' 

CSJ 

A MATRIX 
B MATRIX 

C MATRIX 


D MATRIX 


E MATRIX 


XKT MATRIX 
XMAX MATRIX 

XKS MATRIX 
XLS MATRIX 


ENGINEERING SYMBOLS (REF. 1 ) 
si 
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SUBROUTINE REGMAT 


The segment stiffness matrices, XKS, and the segment load matrices, XLS, are 
passed from SEGMAT to REGMAT via Tapes $2 and and are placed in the XKRTOT 
array and the XLRTOT array, respectively. If kinematic links occur between 
segments in the region, the XKRTOT array and the XLRTOT array are modified 
to represent the situation. In the case of discrete rings the routine RINGER 
is called and provides the necessary matrices. 

A horizontal and vertical partitioning of the XKRTOT array occurs while the 
XLRTOT array is subjected to a horizontal partitioning only. Appropriate 
matrix operations are performed upon the partitions of each array, thus 
reducing the size of the region stiffness and load matrices and resulting 
in increased program capacity. The results of these manipulations are the 
region stiffness matrix, XKR, and the region load matrix, XLR. In passes 
other than the first pass, first cycle, the load calculations are not performed. 

Subroutines Called from REGMAT 

Subroutine SWITCH : Is a routine used to arrange a matrix in a form convenient 

for use by another routine employing a positive definite method for solving 
linear algebraic equations. 

Subroutine CHASE : Is a routine used to obtain the solution X of the linear 

system AX = Y, given at least one right side of Y and the positive, definite, 
symmetric, real coefficient matrix A. 

Subroutine FUTILE: Is a routine called from CHASE and used to obtain the 

factorization of the positive definite, real, symmetric matrix A into the 
product of a lower triangular matrix and its transpose by utilizing a 
Cholesky decomposition. 

Subroutine TRIEQ : Is a routine called by CHASE or EIGYAL to solve a triangular 

system of algebraic equations. 


i 
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FORTRAN CODE 


ENGINEERING SYMBOLS (REF. 1 ) 


SKL MATRIX 
SKLTR MATRIX 

XKRTOT MATRIX 

1 

XLRTOT MATRIX 

SKL22 MATRIX 
REGTOT MATRIX 

STORE MATRIX 

XK11 PARTITION 

XK12 PARTITION 

/ 

X 

XK22 PARTITION 
XK21 PARTITION 


: [ SKL 22] 


" K ll 

i 

1 

k i 2 " 

K 21 _ 

i 

1 

1 

K 22 _ 


L L . 

ft 

N 

! A 

M 

[ 4 ] 
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FORTRAN CODE. 

ENGINEERING SYMBOLS ( REF. l) 


r A 

XL1 PARTITION 

M 


A 

XL2 PARTITION 

H 

A 

r* -l 

XKR MATRIX 

K 1 / 
_ i i 

XLR MATRIX . 

N ■ 
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ENTRY 


REGMAT 
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SUBROUTINE RINGER 


This subroutine reads the discrete ring geometric data, and temperatures, 
and forms the ring stiffness and thermal load matrices. These matrices 
are passed back to either of subroutines REGMAT or STRMAT (see next) as 
necessary, for incorporation into the region or structure matrices, 
respectively. If the present pass is a nonlinear one, the effects of 
prebuckling hoop strain are considered, and the load calculations are 
omitted. 

The calculations in RINGER account for the eccentricity of the ring 
centroid from the base shell wall, and the offset of the ring centroid 
from the shear center. 



FORTRAN CODE 

ENGINEERING SYMOBLS (REF. l) 

HNGTOT MATRIX 

H 

TDEL MATRIX 

H 

KNGLOD MATRIX 

W 

RC 

r 


c 

RS 

r 


s 

XC 

X 


c 

YC 

X 


s 
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FORM RING 
THERMAL LOAD 
MATRIX 


RINGER 
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SUBROUTINE STRMAT 


The region stiffness matrices, XKR, and the region load matrices, XLR, are 
passed from REGMAT to STRMAT via Tape jfb and Tape #8, and are placed in 
the XKSTOT array and the XLSTOT array, respectively. A matrix, BCD, is 
formed to represent the boundary conditions, and, if kinematic links 
occur between regions, the RKL matrix is developed to represent this 
situation. The subroutine RINGER is again called for discrete ring 
matrices. In passes other than the first pass, first cycle, the load 
calculations are not performed. 

As a result of appropriate matrix operations, a reduced structure stiffness 
matrix is formed. In passes other than the first pass, first cycle, STRMAT 
terminates here. In the first pass, first cycle, the following calculations 
are made. Subroutine FLEX, a routine identical to SREVN2 (except for being 
in single precision) with the name changed due to the structure of the 
OVERLAY option, is called to invert the structure stiffness matrix thus 
producing A, the flexibility matrix for the structure. The region end 
deflection array, DRE, is produced as the result of another set of matrix 
operations . 
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FORTRAN CODE 


ENGINEERING SYMBOLS (HEP. 1 


BCD MATRIX 
BCT MATRIX 


[bc] 

[Bcf 


XST MATRIX 


w. 


) 



XKF MATRIX 


A MATRIX 


[ k ]f 

W F 


XSL MATRIX 


N t 


XLS ARRAY 
XFL ARRAY 
DRE ARRAY 


( l )f 

A 

W F 

Hi 
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STRMAT 
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SUBROUTINE INITAL 


As a result of the matrix operations performed in REGMAT, the SKI22, the 
XK2221, and the XK22IP arrays for each region are passed to INITAL. The 
XK1112 and XL1 arrays for each segment, resulting from the matrix procedures 
in SEGMAT, are also passed to INITAL. The region end deflection matrices, 
DRE, which were formed in STRMAT or BCVECT are transmitted to INITAL. 

Following appropriate matrix operations upon these arrays, the force initial 
conditions, the FICS array, and the deflections initial conditions, the 
DICS array, are produced. These arrays combine to form the YICS matrix, 
which contains the true initial conditions for the structure to be analyzed. 

The pertinent counters in the subroutine are: 

NS = segment counter 
NR = region counter 
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FORTRAN CODE 


XK2221 MATRIX 

XK22L2 MATRIX 

DSE ARRAY 
XK1112 MATRIX 

ROTD MATRIX 
DICS ARRAY 
XL1 ARRAY 

ROT MATRIX 
FICS ARRAY 


ENGINEERING SYMBOLS (REF, 1 ) 

AT A 

A , A 

H H 

w 


|s(i)| 

(AD) 

Jift] 
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SUBROUTINE LEBEGE 


The subroutine link LEBEGE receives the YICS array for each segment from 
INITAL via Tape The subroutine FIXEM is called to integrate the 

differential equations of each segment, under true load conditions. 

FIXEM is identical to subroutine SETUP, while WAND corresponds to subroutine 
MAGIC and only consideration of the OVERLAY structure dictates the change in 
names. The subroutines TOBAR, TEMOEG, PLYCO, and PLYNE are similarly 
equivalent to ROBOT, GEOMET, PLICO, and PLINE discussed previously. 

The results of the final integration sequence are the forces and 
deflections at the beginning, intermediate, and end points of each segment. 
In the first pass, first cycle, the prestress and predeforrnation states are 
printed from this routine. In other passes the U, V, W eigenvectors are 
printed, providing convergence has been reached. 
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SUBROUTINES ODE1 AND 0DE2 


Subroutine LEBEGE calls either ODE1 or ODE2 , as necessary, and various geometric 
and trigonometric clues, as well as the predicted values of the variables for the 
differential equations, are passed to this subprogram via label common area LASTEQ. 

The equations in 0DE1 and 0DE2 are identical to those in subroutines DIF1 and 
DIFF2 respectively, with the addition of the auxiliary equations for YAQRH, and 
YAQTH. Subroutines 0DE1 and 0DE2 perform the final integration for each segment 
in the structure utilizing the initial conditions previously obtained, and return 
these values to LEBEGE via label common area LASTEQ. 


The 0DE1, 0DE2 flow charts are identical to the DIF1, DIFF2 flow charts, respectively. 
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SUBROUTINE EIGVAL 

i 


This is the controlling routine for the program eigenvalue calculation 
loop. The necessary matrices for the calculations are passed to COMPAK, 
and thence to routine EIGEN. The controlling routine EIGVAL will 
provide intermediate matrix print as desired. 

Subroutine COMPAK : This routine reads rows of a two-dimensional array 

from a storage unit into core and then stores them into a column Vector. 
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EIGVAL 


ENTRY 


CALL COMPAK: READ IN MATRICES 
& STORE LOWER TRIANGLES IN 
COLUMN VECTORS 


CALL EIGEN 

i 

’ 


RETURN 







oo OOOOOOOO ooooooo ooooooo o o 

O — « ro^^-triNOf^aoa' *-< eg r*i 4* m 'O r» O' o — < cm m o> O 

OO oooooooo - — • ■ — < . — < — < — < • — * . — « r-< r\j r\j r\j r\j r\j f\j eg pn 

oo ooooooo o ooooooo ooooooo oo 

oo oooooooo ooooooo ooooooo oo 

>?■ sr sT'i-»r'T«r<r>r' 4 - <r >r »j- >r >r >r >3- >r >j- ^ >r >r 

m ro foro fr ro m n ^ cn cn m m m f*i r* re ro ro pn 





oooooooooooo 

lAOMDO'OH 

oooooooooo— 

oooooooooooo 

oooooooooooo 

rororor^r^rof^r^. 




SUBROUTINE EIGEN 


EIGEN is the main eigenvalue and eigenvector calculation routine. 

Printout of the eigenvalues and appropriate eigenvectors is provided by 
this routine. First EIGEN prepares the matrices for eigenvalue extraction, 
calling FUTILE, DAGGER, and SWITCH. Then the matrices are passed to 
SYMOEIG, which obtains the eigenvalues and specified vectors using the 
Householder technique. 

Subroutine DAGGER : This routine performs the reduction of the two matrix 

eigenvalue problem to a standard eigenvalue problem. 

Subroutine SYMEIG (Alternate entry point SYMVEC) : This routine obtains 

the eigenvalues and vectors of a real, symmetric, matrix. 

Subroutine TFORM: This routine performs the Householder similarity 

reduction of a real, symmetric matrix to tri-diagonal form. 

Subroutine STURM: Performs the Sturm-sequence calculations of eigenvalues 

of a real, symmetric, tri-diagonal matrix. 

Subroutine PREP (alternate entry point PET) : Determines the number of roots 
greater or equal to a stipulated value. 

Subroutine QSVFC (alternate entry point QWIEL) : Obtains eigenvectors of a 

tri-diagonal matrix and back transforms them to eigenvectors of the original 
real, symmetric matrix. 

Subroutine RANDOM : Random number generator. This subroutine is an MSFC 

system subroutine and the listing is not included here. Most any random 
number generator may be substituted or the argument may be set to zero. 

Subroutine ANDD: This routine "ands" or "ors" the corresponding bits of 

two given variables with each other. 

Subroutine DOTPRO : Performs inner product accumulation using labeled 

common area INFO. 
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SUBROUTINE DETERM 


This is the controlling routine for the program determinant evaluation 
loop. The matrix in question is passed to DET2 where the determinant 
evaluation calculations occur. 

Subroutine DET2 : This routine computes the determinant of a square matrix. 

The actual determinant value is normalized to ±1. 

Subroutine DC0MP2: This routine performs the triangular factorization of 

a square matrix. 

Subroutine SUPER: This routine performs inner product accumulation. 



DETERM 
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DETERMINANT 0F A SPARSE IN-C0RE REAL MATRIX, WITH PR0TECT I 0N 
AGAINST 0VERFL0W AND UNDERFL0W. DETERMINANT IS GIVEN IN A 
SPECIAL F0RM T0 ALL0W REPRESENTAT I 0N 0F VERY LARGE 0R SMALL 
VALUES. 
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SUBROUTINE BCVECT 


This routine converts the eigenvector obtained during a converged 
eigenvalue pass into the form accepted by INITAL. The calculations 
are identical to those in the latter part of STRMAT, and result in 
the DRE array. 



BCVECT 
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SUBROUTINE ETRAP 


This is an error trap routine which can be called by the MAIN routine 
at various stages of program execution. If the indicator NIX is not 
equal to zero, MAIN will call ETRAP and indicate the proper message 
to be printed. 
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